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\0 ■ ABSTRACT 

Cn I The majority of Globular Clusters show chemical inhomogeneities in the composition 

■ of their stars, apparently due to a second stellar generation in which the forming gas 

(— I ' is enriched by hot-CNO cycled material processed in stars belonging to a first stellar 

Oj. generation. Clearly this evidence prompts questions on the modalities of formation of 

I ' Globular Clusters. An important preliminary input to any model for the formation 

2 I of multiple generations is to determine which is today the relative number fraction 

-*— > . of "normal" and anomalous stars in each cluster. As it is very difficult to gather 

5G ' very large spectroscopic samples of Globular Cluster stars to achieve this result with 

i_i| good statistical significance, we propose to use the horizontal branch. We assume 

that, whichever the progenitors of the second generation, the anomalies also include 
^~^ ' enhanced helium abundance. In fact, helium variations have been recently recognized 

^ , to be able to explain several puzzling peculiarities (gaps, RR Lyr periods and period 

distribution, ratio of blue to red stars, blue tails) in horizontal branches. We summarize 
previous results and extend the analysis in order to infer the percentage in number of 
^j^ . the first and second generation in as many clusters as possible. We show that, with few 

• ' exceptions, approximately 50% or more of the stars belong to the second generation. 

tj^ I In other cases, in which at first sight one would think of a simple stellar population, we 

give arguments and suggest that the stars might all belong to the second generation. 

We provide in Appendix a detailed discussion and new fits of the optical and UV data 

of NGC 2808, the classic example of a multiple helium populations cluster, consistently 

^ ' including a reproduction of the main sequence splittings and an examination of the 

l> I problem of "blue hook" stars. We also show a detailed fit of the totally blue HB of 

'^ , M 13, one among the clusters that are possibly fully made up by second generation 

j^ ■ stars. We conclude that the formation of the second generation is a crucial event in the 

life of globular clusters. The problem of the initial mass function required to achieve 
the observed high fraction of second generation stars can be solved only if the initial 
cluster was much more massive than the present one and most of the first generation 
low mass stars have been preferentially lost. As shown by D'Ercole et al. by modelling 
the formation and dynamical evolution of the second generation, the mass loss due 
to the explosions of the type II supernovae of the first generation may be the process 
responsible for triggering the expansion of the cluster, the stripping of its outer layers 
and the loss of most of the first generation low-mass stars. 
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1 INTRODUCTION 

The observations of Globular Cluster (GC) stars are still to 

* E-mail: dantona@oa-roma.inaf.it (FD); vittoria.caloi@iasf- 1 • ^ . 1 • r n ■ , 1 r at ,11 

^. ,, ,^, be interpreted m a lully consistent frame. Nevertfieless, a 
roma.inaf.it (VC) , . -^ . , r , ^U 

t This work has been supported through PRfN INAF 2005 "Ex- general consensus is emerging on the fact that most GCs 

perimenting stellar nucleosynthesis in clean environments" and can not be considered any longer "simple stellar popula- 

PRIN MIUR 2007 "Multiple Stellar Populations in Globular Clus- tions" (SSP), and that "self-enrichment" is a common fea- 

ters: census, characterization, and origin". ture among GCs. This consensus follows from the well known 
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"chemical anomalies", already noted in the seventies (such 
as the variations found in C and N abundances, the Na- 
O and Mg~Al anticorrelations) . Recently observed to be 
present at the turnofF (TO) and among the subgiants (e.g. 



Gratton et ai] I2001I : iBrilev et al.1 booj |2004| : ICohen et all 
2OO5I ). they must be attributed to some process of "self- 



enrichment" occurring at the first stages of the cluster life, 
as the same authors quoted above suggest. There was a first 
epoch of star formation that gave origin to the "normal" 
(first generation, hereinafter FG) stars, with CNO and other 
abundances similar to Population II field stars of the same 
metallicity. Afterwards, there must have been some other 
epoch of star formation (second generation, hereinafter SG), 
including material heavily processed through the CNO cy- 
cle. This material either comes entirely from the stars be- 
longing to the first stellar generation, or it is a mixture of 
processed gas and pristine matter of the initial star form- 
ing cloud. We can derive this conclusion as a consequence of 
the fact that there is no appreciable difference in the abun- 
dance of elements such as Ca and the heavier ones between 
"normal" and chemically anomalous stars belonging to the 
same GC. Needless to say, this statement does not hold for 
u Cen, which must indeed be considered a small galaxy and 
not a typical GC. In the following, we will only examine 
"normal clusters" , those which do not show signs of metal 
enrichment due to supernova ejecta. The homogeneity in the 
heavy elements is an important fact that tells us, e.g., that 
it is highly improbable that the chemical anomalies are due 
to mixing of stars born in two different clouds, as there is 
no reason why the two clouds should have a unique metal- 
licity. In addition, the clusters showing chemical anomalies 
have a large variety in metallicities, making the suggestion 
of mixing of two different clouds even more improbable. The 
matter must have been processed through the hot CNO cy- 
cle, and not, or only marginally, through the helium burning 
phases, since the sum of CNO elements is the same in the 
"normal" and in t he anomalous stars (e.g. Smith et al.ll 19961: 
llvans et all 1 19991 : ICohen fc Melendej 120051 ) . ICarretta et al.1 
((20051) find that actually the CNO is somewhat -but not 
much- larger in the SG stars of some GCs. Therefore, the 
progenitors may be either massive asympt otic giant branch 
(AGB) stars ( e.g. I Ventura et al.ll20oil . i2002l flor fast rotating 
massive stars IjDecressin et al.ll2007^ . In both cases, models 
show that the ejected material must be enriched in helium 
with respect to the pristine one. The higher helium con- 
tent has been recognized to have a strong effect on hori- 
zontal branch (HB) morphology, possibly helping to explain 
some features (gaps, hot blue tails, se cond parameter) which , 
until now, have defied explanation (JD'Antona et al.ll2002h . 
Along these lines, a variety of problems has been examined: 
the extreme peculiarity of the HB morphology in the mas- 
sive cluster NGC 2808, l|D'Antona fc Caloill2004l): the sec- 
ond p arameter effect in M 13 and M 3 ( Caloi fc D'Antonal 
[2005); the peculiar features in t he RR Lyr v ariables and 
HB of NGC 6441 and NGC 6388 (|Caloi fc D'Antonaii2007ii) . 



The presence of strongly enhanced helium in peculiar HB 
stars has been confirmed, f or NGC 2808 and NGC 6441, by 
spectroscopic observations (jMoehler et al.ll2004l : lBusso et al.l 
[2007,' ). 

Beside this spectroscopic evidence, an unexpected fea- 
ture has recently appeared from photometric data: the 
splitting of the main sequence in NGC 2808. After first 
indic ations from a wider t han expected colour distribu- 
tion (JD'Antona et al.l |2005| ). recent HST observations by 
IPiotto et al.l ( 20071 ) leave no doubt that there are at least 
three different populations in this cluster. This came af- 
ter the first discovery of a peculiar blue main sequence in 
LO Cen (JBedin et al.ll2004l ), interpreted again in terms of a 
very high helium content (|Norrisl [JOOJ : IPiotto et al.ll2005l ). 
The above mentioned cases can be considered as "extreme" 
ones, in the sense that no explanation had been attempted 
for them before the hypothesis of helium-enriched popula- 
tions. Less critical situations, such as the HB bimodality in 
NGC 1851 and 6229, had been tentatively explained in terms 
of a unimodal mass distribution with a large mass disper- 
sion (0.055 - 0.10 Mq, Catelan et al. 1998). But even such 
a rather artificial assumption could not help in the case of 
NGC 2808, which, as hinted before, finds a possible solution 
only in terms of varying helium in multiple stellar gener- 
ations. Therefore, we consider appropriate to apply to the 
less peculiar cases the solution found plausible for the most 
peculiar ones. In fact, notice that split main sequences and 
strong bimodalities are only the tip of the iceberg of the 
self-enrichment phenomenon. In most clusters the higher 
helium abundances remain confined below Y~ 0.30, and 
the presence of such stars would not be put in evidence 
eithe r from main seque nce observations (JD'Antona et al.l 
|2002| : ISalaris et al.l I2OO6I ) , or from a naif interpretation of 
stellar counts on the HB, as we shall discuss in Sect. 14.11 

If we wish to shed light on the entire process of forma- 
tion of GCs we must have a rough idea of the total number 
of SG stars. We will then analyze the HB in terms of pop- 
ulations differing in Y, using one or more of the following 
peculiar features: 

(i) bimodal HBs and HB gaps; 

(ii) presence of blue HB stars and very long period RR 
Lyr's in high Z clusters; 

(iii) peaked number vs. period distribution of RR Lyr's; 
(iv) blue-HB clusters. 

In this paper we show the results of such an interpretation 
for several GCs. We start from a reanal yis of the NGC 2808 
data, taking into account the results bv lPiotto et al.l (|2007l ) 
for the main sequence , and the ultraviolet HST data by 
ICastellani et al.l (|2006l ): we summarize the results already 
published and discuss briefly the other clusters. The table 
of the derived FG and SG percentages is the basis to discuss 
the clusters' dynamical evolution required to produce the 
high fraction of stars presently belonging to the SG. 



1 If tlie ICarretta et al.l ll2005h CNO data really indicate that a 
limited number of third dredge up (e.g. Ilben &: Renzinil 119831) 
episodes plays a (small) role in the nuclear processing of the mat- 
ter giving origin to the SG, the massive AGB progenitors are 
possibly favoured. 



2 THE BASIC MODEL 

Here we summarize why and how a helium enri chment mod- 
ifies the HB morphology (JD'Antona et al.ll2002l ) , and the ba- 
sic inputs of the HB and main sequence (MS) simulations 
adopted to constrain the FG and SG. 



2.1 Red giant mass, mass loss, and helium content 

The evolving mass in a GC is a function of age, metal- 
licity and helium content. From the well defined turnoff 
of GCs, it is evident that any age spread must be much 
smaller than the global age (10-13 Gyr). In addition, the 
metallicity spread among cluster members are contained 
within ~ 0.04dex Q Thus, once fixed age and metallicity, 
the evolving mass is a function of the helium content only, 
and it decreases when helium is increased. The dependence 
5Mrg/5Y ~ — I.3M0 is such that, e.g. the small increase 
in helium content from the primordial value Y=0.24 to the 
moderately higher Y=0.28 decreases the evolving mass by 
~O.O5M0. 

During the latest phases of red giant evolution, both 
normal and helium enhanced stars lose mass. We assume 
that mass loss follows Reimers' law (|Reimerall975l ') 

M« = 4.10-i^«^ (1) 

where rjn is a free parameter directly connected with the 
mass loss rate and L, R and M are luminosity, mass and ra- 
dius expressed in solar units. This expression has no explicit 
dependence on helium abundance or metallicity. While the 
independence from the metallicity may be questioned, the 
helium content, at the level of variation we are considering, 
should not affect the mass loss rate, as helium has no strong 
eflfect on surface opacities and possible grain formation. It 
turns out that also the total mass lost by giants with differ- 
ent hefium content and similar age does not depend on the 
helium content. In fact, we computed tracks of stars with 
different helium content, and having different mass so that 
they have the same evolving age, and we find that the total 
mass lost differs only by 0.001 - 0.002 Af0when the models 
reach the helium flash. On the other hand, the helium core 
mass at the helium flash depends on the helium content, 
decreasing with increasing helium, but at a much smaller 
rate than the decrease of the red giant mass. Therefore, the 
ratio of core mass to the remnant mass is larger for larger 
Y, and these stars will occupy a position on the ZAHB at a 
bluer colour than stars with lower Y. This means that, if the 
cluster contains stars of FG with "standard" Y, and stars of 
SG with larger Y, these latter will be "bluer" than the FG 
stars. 

Assuming the point of view that the HB contains stars 
with different helium content, the estimate of age and metal- 
licity on the one side, and the HB morphology on the other, 
will indicate which part of the HB population belongs to 
the FG. This component should have a uniform helium con- 
tent, probably c lose to the Big Bang abundance (e.g. Y~0.24 
ICoc et aLll2004l ). For the FG of the most metal rich clusters 
we will assume a larger initial, uniform helium Y=0.25. On 
the contrary, the SG will most probably show a spread in Y, 
for two main reasons: the self-enriched material in fact 1) 



may come from different progenitors, having different chem- 
ical peculiarities or 2) it may be diluted in different fractions 
with matter from the FG. In both cases the helium abun- 
dance may differ among the SG stars. Notice however that if 
we have physical reasons (e.g. based on one of the peculiari- 
ties listed in Sect. QJ) to attribute to a star a helium content 
larger than the helium of the FG, we know that this star he- 
longs to a SG, even if its Y is not much larger. Of course, 
the total amount of SG self-enriched gas differs if the de- 
rived Y does, or does not, result from dilution with pristine 
matter. 



2.2 The grids of HB models 

The basis of the synthetic HB distributions are stellar 
models computed w ith t he code ATON2 . 0, de scribed in 
IVentura et all ([Hii) and iMazzitelli et all (| 19991) . The HB 
models have been evolved until the disappearance of he- 
lium in the convective core. We have adopted metallicity 
Z=2 10"* for the metal poor GCs; Z^IO"^' and Z=2xl0"^ 
for the intermediate metallicity clusters, which represent the 
majority of data sets, and Z=6xl0~^ for the high metallic- 
ity clusters. The helium core mass of the models is set at 
the helium flash core mass of the previous evolution, de- 
termined by evolving models for each couple (Z, Y) for an 
age of 11 X 10^ yr. The metallicities of individual clusters may 
be slightly different from those adopted here, but the main 
aim — to distinguish between normal-helium and enhanced- 
helium stars in the construction of the HB — can be satis- 
factorily achieved. We will see that the ratios FG/SG are 
well defined from gross characteristics of the HB, and not 
by minute details. We computed grids of models for Y=0.24, 
0.28, 0.32 and 0.40. 

We computed HB models up to Toff~31000K, that is 
the usually accepted limit for standard ZAHB models re- 
sulting from a He-flash occurring at the red giant branch 
tip. 

These models can not explain the extreme "blue 
hook" stars (Teg up t o ~37000K) present in o^Cen M54 
l|Rosenberg et al.ll2004h, NGC6 388 ||Busso et al.ll2007l ') and 
NGC 2808 (|Moehler et al.ll2004l '). and generally explained as 
a result of the mixing of processed matter with the very 
small residual hydrogen envelope, consequent to a late ig- 
nition of the helium flash, along the white dwarf cooling 
sequence l|Sweigartlll997l : iBrown et al.ll2001^ . Mixing raises 
the helium (and ca rbon) abundance in the envelope (see also 
ICassisi et a l. 2003) and the st a r sett les at larger Teff. In 
addition, ID'Antona fc Ventural l|2007l ) have suggested that 
deep mixing occurs, independently from a late helium flash, 
along the RGB evolution of giants belonging to the very high 
helium population (Y~0.35 - 0.40) as deep mixing in this 
case is not forbidden by the molecular weight discontinuity 
barrier. So, deep mixing may involve most of the very he- 
lium rich stars, and increase, even considerably, their surface 
helium abundancqj. 



^ According to lGratton et al.l 1120041 ) : "Low upper limits in the 
spread of the abundances of Fe (of the order of 0.04 dex, r.m.s.) 
have been found for several clusters from both spectroscopy and 
from the widths of main sequence (MS) and red giant branch 
(RGB) stars in the colour magnitude diagrams CMDs (for sum- 
mary and discussion, see Suntzeff 1993). At present the verdict 
on Fe variations in CGs except u) Cen must remain 'not proven'." 



^ The hypothesis of deep mixing is useful to increase the proba- 
bility of obtaining extreme HB and blue hook stars. If we use only 
the very late flash hypothesis, just a few stars could happen to 
have the appropriate envelope mass, such that they do not ignite 
the flash on the RGB, or do not leave a He-white dwarf remnant. 
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In order to simulate the possible late-flash induced mix- 
ing (and any other kind of deep mixing), for the set of 
Z^O.OOl adopted to simulate the HB of NGC 2808, we com- 
puted tracks of very low masses for Y=0.45, 0.50, 0.60, 0.70 
and 0.80, having core mass M=0.4676 Mq, the helium core 
mass at the flash for the Y=0.40 models. 



2.3 Main sequence and turnofF models 

For comparison with the main sequ ence colour distribution 
in NGC 2808 bv lPiotto et al.l (J2007l ). we use the isochrones 
for Z=10~^ and ¥=0 .24, 0.28, 0.32 and 0.40 described in 
iD'Antonaet all (|2005l ;). 



2.4 Simulations 

We adopt an appropriate relation between the mass of the 
evolving giant Mrg and the age, as function of helium con- 
tent and metallicity. The mass on the HB is then: 



(120041). The results are very s imilar if we adopt the clas 



van Albada fc Bakerl (|l973l ) relation. The real problem 



Mhb = Mrg{Y,Z)~/\M 



(2) 



AM is the mass lost during the RG phase. We assume that 
AM has a gaussian dispersion o around an average value 
AMo and that both AMo and a are parameters to be de- 
termined and do not depend on Yy Therefore the HB mass 
varies both due to the mass dispersion g around the average 
assumed mass loss, and due to the dependence of the RG 
mass on the helium content. As the evolving mass decreases 
with increasing helium content, the stars with higher helium 
will populate bluer regions of the HB. In each cluster, we 
must identify the FG section of the HB (e.g., in NGC 2808 
this is easily identified with the red clump). Then we assume 
a primordial (Y=0.24) or a minimum (Y=0.25) helium con- 
tent for the fraction of HB stars considered to belong to the 
FG. This fraction can be adjusted in order to reproduce the 
features we attribute to the FG (e.g., the peaked RR Lyr 
period distribution in M3, see later). We then assume that 
the rest of HB stars has larger Y, and adopt different N(Y) 
distributions in order to reproduce the other parts of the 
HB number vs. colour distributions, when available. 

We fix AMo and a, and extract random both the 
mass loss and the HB age in the interval from 10 yr to 
lO^yr, according to the chosen Y distribution. We thus lo- 
cate the luminosity and Tea along the evolution of the HB 
mass obtained. We identify the variable stars as belong- 
ing to a fixed TeS interval and comp ute their period ac- 
cording to the pulsation equation (1) bv lDi Criscienzo et al.l 

* A complication in the interpretation of HB morphologies, not 
included in this work, arises if not only the helium content be- 
tween FG and SG varies, but also the total CNO content. This 
additional problem comes out from the observations of NGC 1851 
(see Sect. 3.2 ) whose HR diagra m show s a splitting in the sub- 
giant branch JMilone et al.ll2008l ). While ICassisi et al.l ll2008i') are 
able to explain this feature by assuming an SG with about double 
CNO content, and same age as the FG (thus in the classical self- 
enrichment sce nario), the subseque nt interpretation of the HB 
morphology by ISalaris et al.l 1120081) can not be simply achieved 
by doubling the CNO, but it requires some extra mass loss for 
the SG stars. Here, our main aim is to find out the percentages of 
FG and SG, and this can be obtained also within the framework 
of our simpler assumptions. 



is given by the choice of the exact boundaries of the RR 
Lyr strip, that affect strongly the number and mean pe- 
riod of the RR Lyrae v ariables (see, e.g., the discussion in 
ICaloi fc D'Antonall2008l l. 

The L and Tefi values are transformed into the differ- 
ent observational bands. As most observations are available 
in the B and V bands, and/or in the Bessell's I, we derive 
the visual magnitude Mv and the B-V or V-I colours, us - 
ing the transformations bv lBessell. Castelli. fc Plea ([1993). 
Although our main aim — to understand the different frac- 
tions of FG and SG stars — is not affected by this problem, 
we remark that the B-V and V-I colours saturate at large 
Tcff. The bolometric corrections become very large, so the 
number vs. magnitude distribution may suffer some uncer- 
tainties, which can be avoided by using different magnitudes, 
e.g. the ultraviolet HST bands. We exemplify such compar- 
isons for the case of the clusters NGC 2808 and M 13 in 
the Appendix A. T he transformations for the ACS - HST 
bands are taken from'Bedin et al.' ('2005'). The WFPC2-HST 
relations are by[Origlia & Leithcrcr (2000), plus additional 
transformations kindly provided by L. Origlia. 

We associate a gaussian spread in colour and magnitude 
to each point, in order to simulate the impact of observa- 
tional errors. We do not include binaries in the simulations. 
We mostly compare the theoretical simulations and the ob- 
servations by looking at the number counts vs. colour in the 
horizontal part of the HB (e.g. for the red part and the RR 
Lyr) or vs. magnitude, for the vertical part, if the blue HB 
is very extended. 



3 BIMODALITY AND GAPS IN THE HB 

There is a huge literature which has defined and attemped 
to u nderstand the gaps on the blue side o f the HB (see, 
e.g. iFerraro et al.l 119981 : iPiotto et al.l ll993. and references 
therein). While a gap at Teff~ 10*K sh ould probably be at- 
tributed to the operation o f diffusion (JGlasoev et al.lll989l : 
ICrundahl et al.l 19991 : ICaloilll999i ). other gap s may have to do 
with discontinuities i n the helium content (JD'Antona et al-l 
l2005l : lLee et al. ' 2005). In particular, a bimodal distribution 
in colour characterizes the HB of some clusters, that have 
well populated blue and red sides of the HB, with limited 
or null population of R R Lyr variables . NGC 2808 is a pro- 
totype of this class, and lCatelan et al.l ([1993) were not able 
to interpret it in terms of a unique mass distribution, even 
with a mass spread as large as 0.3 M©. A multimodal mass 
distribution was then required to explain this HB. 

The HB bi modality in NGC 2808 was the main hint 
used by iD'Anto na & Caloi (2004) to infer the presence of 
multiple stellar generations differing in Y in the cluster, an 
interpretation nicely supported by t he subsequent observa- 
tion of the main sequence splitting ([D'Antona et al.ll2005l : 
IPiotto et al.ll2007r ). In addition, in this cluster are found the 
still misterious "blue hook" stars. So NGC 2808 appears 
as an ideal benchmark for the application of the multiple 
population hypothesis. We summarize in the following our 
present understanding of its modeling, and extend it to other 
situations. 



3.1 NGC 2808 

If age, metallicity and mass loss are such that normal- 
hehum, FG stars, populate a red clump, the SG stars with 
helium enhancement will tend to populate the bluer HB and 
the RR Lyr region. If there is a gap between the normal- 
helium stars and the minimum helium content of the second 
generation, the case of NGC 2808 may appear: a red clump 
(FG), almost no RR Lyr (due to the helium gap) and a blue 
HB with la rger helium cont ent (star ting from Y^ 0.28 ac - 
cording to iD'Antona fc Calo i 2004: iD'Antona et al.ll2005h . 
In those papers, we had modelled the mass loss along the 
RG branch by assuming that the larger Y would provoke 
a slightly larger global mass loss, as the evolving giants 
with higher Y are loss mass ive, and t hus have smaller grav- 
ity (see, e.g. iLee. Demargue. fc Zinnlll994 ). We have ascer- 
tained that this is not the case (Sect. I2.1|l . so the HB has 
to be modeled by assuming the same average mass loss for 
both normal Y and higher Y red giants. In addition, we 
try to model better t he EBT2 and EBT3 (in the definition 
by iBedin et al.l 2004h blue clumps, which contain the ex- 
treme HB (JD'Cruz et aTlbOOd: iBrown et al.l^2001^ and the 
"blue hook" stars ( Moehler et al.l 120041 ') . respectively, and 
compare also simulations based on HST ultraviolet and vi- 
sual bands. The details can be found in the Appendix. We 
show that th e new simulations a re consistent with both the 
triple MS bv lPiotto et al.1 (120071 ') and the HB star distribu- 
tion. The main differe nce with respect to the analysis by 
ID'Antona et al.l l|2005i ) is that the intermediate Y popula- 
tion is now clustered at Y~0.31. The cluster again results 
divided into 50% normal-helium stars, and 50% helium en- 
riched stars, although the very high helium (Y~0.385) stars 
are only ~15%. 



3.2 NGC 1851, NGC 6229 

Following the guidelines mentioned in Sect. [T] we take the 
point of view of interpreting clusters with a bimodal HB 
in terms of multi ple populations , as in the "pa radigmatic" 
case o f NGC 2808. ICatelan et al.l (Il998l ') and lBorissova et al. 
l|l999h describe the complex HB structure of the cluster 
NGC 1851 and NGC 6229, respectively: bimodal, with few 
RR Lyr variables, a gap on the blue HB and possibly some 
extreme blue HB members, at the luminosity of the turn- 
off. In NGC 6229, the number ratio of the red, variable and 
blue HB components are: B:V:R=0. 59:0. 08:0. 33. The RR 
Lyr average period is of the Ool type. We assume therefore 
that the RR variables belong to the first generation, since 
they have the period appropriate to the metal and helium 
content expected for these stars. So we expect that roughly 
41% of cluster members belong to the first generation and 
59% to the second, helium-enriched one. The extra- helium 
allows these latter stars to reach the bluer positions beyond 
the variable region, to which they should have been con- 
fined by a chemical composition of Y~0.24 and Z~0.001. 
(See later the case of M3). Of course, not all the helium- 
enhanced stars will be found on the blue, and vice-versa, 
but the distribution will be substantially the one mentioned 
above. 

Similarly, for NG C 1851 we have: B:V:R:=0.30:0 10:0.60 
JCatelan et all Il998l ). or B:V:R~0.32;0.12;0.56 l|WalkeiJ 
119981 ). In this case the RR Lyr are again Oosterhoff type 



I, but their periods are longer, so that the variables may 
belong to the SG. Recently, this latter cluster has been dis- 
cover ed to harbour a double subgiant branch IjMilone et al.l 
I2008I ). that can be interpreted as the presence of an FG 
and SG with differe nt total CNO abundances and same age 
(jCassisi et al.ll20p8|). The "bright" s ubgiant branch contains 
55±5% of stars ( Milone et al.ll2008r ) and may correspond to 



the red part of the HB. Thus the FG should con tain ~55% of 
the to tal cluster stars. As we already remarked, ISalaris et al.l 
(|2008h analyze the HB stellar distribution, and are not able 
to fit the blue and red side of the HB with different helium 
(and CNO-Na abundances) and the same mass loss on the 
RGB. While modelling of this HB may require a database 
of HB tracks computed with different Y and different com- 
positions in CNO, for our present purposes it is sufficient 
to see that the consistent bimodality of both the SGB and 
HB indicate the presence of an FG and SG, with the given 
proportions. 



4 CLUSTERS WITH HIGH METALLICITY 
AND PECULIAR HB 

4.1 NGC 6441 and NGC 6388 

The case of the hig h metallicity cluste r NGC 6441 has been 
fully discussed by ICaloi fc D'Antonal (|2007l ). who showed 
that very helium rich stars are also present among the red 
clump stars. The analysis is able to ex plain not only th e 
anomalous long periods of the RR Lyr (jPritzl et al.ll2003l ). 
but also the extension in magnitude of the red clump (any 
attempt to attribute this thickness to differential reddening 
has failed — [Raimondo et al.l l|2002r )'). and the hot blue side 
of the HB. Therefore, the morphology requires helium en- 
richment not only for the bluer side of the HB, as we could 
naively think, but even for the red clump stars. The physical 
reasons for this interpretation is the following: helium core 
burning stars havi ng high Y and Z make lo ng loops from red 
to blue in the HB l|Sweigaft fc Gro ss 1976). In fact, both the 
higher mean molecular weight -leading to a high H-burning 
shell temperature- and the high metallicity -leading to a 
stronger CNO shell- conspire towards the result that the 
H-shell energy source prevails with respect to the He-core 
burning. The consequent growth of the helium core leads the 
evolution towards the blue. Therefore, if we must explain the 
luminous (long period) RR Lyr by stars having high helium, 
the same stars will also populate the red clump: this is ex- 
actly what we observe: if the helium content is not as large 
as Y~ 0.35 in the red clump, the HB finds no satisfactory 
explanation. The percentage of helium enriched stars is in 
this case ~60% for NGC 6441. 

We performed a similar analysis for NGC 6388 (see 
Fig. 1). The main difference among the two clusters is that 
NGC 6388 seems to have a higher tail of very high helium 
(Y>0.35) stars, reaching ~20%. 

The re are other analys es for these clusters in the lit- 
erature: Buss o et al.l (|2007l ) consider as peculiar only the 
blue HB stars. T his would limit the SG to ~ 15%. Also 
lYoon et al.l (|2008h attribute the presence of a SG to the RR 
Lyr and hotter stars only. They seem to be able to obtain 
the RR Lyr long periods with only Y~0.3, but we have no 
details about their models to understand this difference. Our 
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Figure 1. Synthetic HB simulation for NGC 6388. A similar 
analysis for NGC 6441 is shown in Caloi & D'Antona (2007). The 
HB data are taken from Piotto et al. (2002). In black we have the 
observations, while red and blue are the simulated stars. In red are 
the Y=0.25 stars, in blue those with Y>0.25 for a total of 1300 
stars. The observed clump distribution (full line histogram) is 
shown on the right, superimposed to the theoretical distribution 
(dash-dotted histogram). The bottom panel shows the number 
vs. helium distributions assumed for the simulation. The number 
of stars with primordial helium Y=0.25, N(0.25), is indicated in 
the label. 



models require Y up to ~0.35 to fit the long periods of the 
RR Lyr, but, as explained above, this high Y also helps to 
reproduce a thickness > 0.7 mag of the red clump. On the 
other hand. lYoon et al.l (|2008 ) attribute the thickness of the 
red clump to differential r eddening, a hypo t hesis in contrast 
with the data analysis bv lRaimondo et al-l l|2002l ). 



4.2 47 Tuc 

In 47 Tuc, the prototype of metallic GCs, the red clump 
is much less thick in magnitude than in the two anomalous 
clusters above. In Fig. 2 we show the histograms of the num- 
ber of stars in the red clump of NGC 6441, NGC 6388 and 
47 Tuc versus magnitude. The magnitudes have been nor- 
malized so that the peak in the distribution coincides for 
all the clusters. The "thickness" in magnitude of NGC 6388 
and NGC 6441 is larger than for 47 Tuc. The excess of stars 
at smaller luminosities below the maximum is probably due 
to the larger observational errors, but the (asymmetric) ex- 
cess at higher luminosities is most easily interpreted as due 
to stars with helium much higher than normal. Based on 
this feature only, we infer that the 47 Tuc SG should not be 
larger than ~25% of the stars. Observations show that CN 
strong and CN weak stars in the cluster are about in similar 
percentages ( |Brilev et al.ii2004 ). and if these two groups are 
to be interpreted as FG and SG, we have a contradiction. 
An escape from this problem can be found if the first stel- 
lar generation in 47 T uc has a larger initial helium content 
l|Salaris fc Weisslll998l ). 



Figure 2. Plot of the observed number of stars versus magnitude 
for the red horizontal branch of three metal rich clusters: NGC 
6441 (dots), NGC 6388 (full line), and 47 Tuc (dashed line). See 
text for the interpetation. 



5 THE PEAKED PERIOD DISTRIBUTION OF 
RR LYR STARS 

The prototype of HBs has often been considered the HB of 
the cluster MS: it is well populated both in the red part, in 
the RR Lyr's and in the blue side, without a blue tail. The 
HB distribution among red, variable and blue members can 
be reproduced by assuming an average mass loss along the 
RGB, with a standard deviation a ^ 0.025 Mq. 

However, there are two important facts to be men- 
tioned: i) the detailed colour distribution along the HB is 
by no means uniform, and ii) the RR Lyr period distribu- 
tion appears strongly peaked, a feature that cannot be un- 
derstood in te rms of a more or less uniform mass distribu- 
tion (C astclla ni fc Tornambd |l98ll: iRood fc Crocked Il989l : 



ICatelanll2004l: ICastellani et al 



2001 



We have examined in 
since there are 



detail this case l|Caloi fc D'Antona! IJOOJ ): 
several other clusters showing the same problem, we sum- 
marize the main points of the model for M3 and extend the 
interpretation to other clusters. 



5.1 The period distribution of M 3 

As said before, the RR Lyr period distribution in M 3 is 
highly peaked (see Figure 3, left side, full line histogram). 
ICastcllani ot al. (20051) realized that the only way to repro- 
duce this peak was to reduce the dispersion in the mass lost 
in the RGB. Once obtained the correct distribution with a 
given recipe for the mass loss, the authors had to assume 
a different average mass loss, with a different dispersion, to 
account for the blue side of the HB. In the hypothesis of mul- 
tiple helium enhancements, the blue side is naturally popu- 
latcd by helium rich stars, as shown bv lCaloi fc D'Antonal 
(2008). In the simulations, we can explain both the period 
distribution and the colour distribution along the HB for the 
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red, variable and blue regions. The analysis however poses 
another problem: we find that the dispersion in mass loss 
along the RGB must be at most a '^ 0.003 Mq to be con- 
sistent with the period distribution. The question remains 
whether this small dispersion is peculiar to M 3, or rather 
we have always been misled by the overall reproduction of 
the HB morphology, for which previous simulations (at con- 
stant Y) needed a dispersion in mass loss of some hundreths 
of Mq, as quoted before. 



5.2 NGC 3201, NGC 7006, M5 

We examined NGC 3201 according to the same scheme used 
for M 3, o n the basis of th e data by iLavden fc Saraiedinil 
(|2003t l and lPiersimoni et al.1 (J2002l l. Un fortunately, the clus- 
ter is affected by differential reddening (jvon Braun fc Mated 
I2OOII . and references therein). While it is possible to per- 



form reddening estimates for the single RR Lyrae variables 
(Sturch 1966, Blanco 1992; see Layden fc Sarajedini 2003), 
the same is not so easy for non variable stars. In fact, when 
we tried to reconstruct the colour distribution on the HB, we 
found many non variable stars in the variable region. There- 
fore we did not attempt to reproduce the detailed colour 
distribution as we did for M 3, but only the RR Lyr period 
distribution and the overall division among red, variable and 
blue members. 

We compared the period distributions in M 3 and in 
NGC 3201, using both the data by Piersimoni et al. and 
Layden fc Sarajedini. Notwithstanding the large difference 
in the total numbers of variables with an established period 
(more than 200 variables in M 3, slightly more than 50 in 
NGC 3201), one finds a strong similarity in the period dis- 
tributions (see Fig. O left side). So we have again a peaked 
distribution, and even a dip at the same period! Therefore 
the NGC 3201 solution is quite similar to the one for M 3. 
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Figure 3. In the left figure we plot the number of RR Lyr vs. pulsation period for NGC 3201 from two different data bases: dashed line, 
by L ayden et al. 2003, and dotted by Piersimoni et al. 2002. We also plot the period distribution for M 3 (full line ICorwin fc CarnevI 
I2OOM . but these data have been divided by a factor two, in order to allow an easy comparison with NGC 3201. In spite of the different 
total numbers, the distributions are very similar. On the right we plot the NGC 3201 period distribution by Piersimoni et al. (dotted) 
and its simulation (full line), discussed in the text. 
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Figure 4. Period distribution fits for M 68 (left) and M 15 (right). Dotted histograms show the observed data, full line histograms are 
the simulations. 



Our best simulations give 37% of HB members with helium 
enhanced from Y=0.26 to 0.28, that is, coincident with the 
percentage of the blue HB stars. The mass dispersion op- 
timum is even smaller than in M3 {a ~ 0.0015 Mq). The 
period simulation is shown (full line histogram) in the right 
side of Fig. [3] 

What said for NGC 3201 can be repeated for NGC 7006, 
whose RR Lyrs sh ow again a peaked period distribution 
IjWehlau et al.ll 19991 ) . Given the morphology of the HB, com- 



posed mainly by red and variable stars, the required percent 
of helium enhanced objects is of about 27%, up to Y about 
0.27. 

The case of M 5 (NGC 5904) appears less simple, since 
the RR Lyrs period distribution presents two peaks. We 
can only approximate this distribution, while we succeed in 
reprodu cing the detailed co l our d istribution along the HB 
given bv lSandauist fc Boltd l)2004h . Within these limits, we 
estimate 70% of SG stars, up to a helium content of about 



0.31 (we remind that M5 has more blue HB members than 
M 3, with the peak of population in the blue). 



5.3 Very metal poor clusters: M 15 and M 68 

A very good fit of the RR Lyr peaked period distribu- 
tion in M 15 and M 68 is obtained with the same method 
adopted for M 3, but obviously using the tracks of metal- 
licity Z=2xl0~^, adequate to describe these two clusters. 
The resulting period distributions are in Fig. 2] We must 
say that the fit of HR diagram, on the contrary, is much less 
successful than for the other clusters we have analyzed so 
far. The discrepancy is the following: the blue part of the 
HB in the simulations is achieved with Y in the range 0.26- 
0.28 for M 68, and 0.26-0.30 for M 15, but the luminosity of 
the blue side is too large with respect to the data. We sug- 
gest that helium is not the only parameter varying in these 
clusters, and that also the total CNO abundance should be 
varied, giving origin to a more complex description, still to 
be explored (see, e.g., the problem of NGC 1851 described 
above) . 



6 THE "BLUE-HB" CLUSTERS 

A specially intriguing case is presented by the cluste rs with a 
preva lently blue HB, that is, with a HB type 1-2 (|Dickena 
119721 '). Almost - if not all - of them fall in the category 
of the clcissical "second parameter" problem. A fact that is 
not always considered as it deserves is that these clusters 
represent the ma jority of the populat ion of the intermediate 
metallicity GCs (|Alcaino et al.lll999r ). So these clusters are 
not an exception, but rather the rule, that is, the most com- 
mon result of the GC formation process. We shall examine 
some of these cases. 



6.1 M 13: multiple populations, but all belonging 
the SG? 

We analysed in detail the case of M 13, finding substantial 
support to the hypothesis that only the second star genera- 
tion has survived in the cluster. This conclusion was reached 
on the basis of the relative positions of t he turn-off, the red 
giant bump and the horizontal branch (|Caloi fc D'Antonal 
I2OO5I ). This result, if confirmed, would be the first direct ev- 
idence for the presence of a substantial helium enhancement 
in GC stars. In the Appendix B we present a detailed analy- 
sis of the H B of this clust e r, in t he HST plane F555 vs. F336- 
F555 from lFerraro et al.l (|l998l '). The helium distribution is 
compared with the one necessary to reproduce the blue HB 
of NGC 2808, with the same assumption on age, metallic- 
ity and mass loss on the RGB. The comparison shows that 
M 13 stars have a shallower distribution in helium, but may 
have a peak at Y~0.28, and a similar peak at Y~0.38. In 
addition, it com pletely lacks th e red c lump stars, that is the 
FG stars. In the lFerraro et al.l (|l998) sample we use, there 
is a total of 221 HB stars, populating the upper, medium 
and low luminosity blue HB. According to the chosen simu- 
lation, a fraction 40% of stars has Y=0. 28-0.29, another 30% 
has Y~0.31-0.35, while the low HB is reproduced by taking 
~30% of stars at Y=0.38. Thus the group of stars mostly 



contribut ing to the RG bump is tha t at Y=0.28-0.29, as sug- 
gested in lCaloi fc D'Antonal l|2005l ). The conclusion is that 
in M 13 there are multiple populations, but, according to 
this interpretation of the data, they all belong to the SG! 

On the other hand, several observations show the exis- 
tence of chemically normal M 13 mem bers, for example in 
Na abundance (e .g. J ilacho wski et al.l ll996'l and in C and N 
abundances (e.g. iBrilev et al.l '2004l. However, if the Y~0.28 
population is the result of the dilution of pristine cluster 
matter (Y=0.24) with highly Y enriched ma tter (Y > 0.34) 
as envisa ged in many formation scenarios i Decrcssi n et al.l 
!2007; Ven tura fc D'Antonall2008bl : lD'Ercoie^ ct al. 200^), not 
necessarily the main chemical anomaly indicators will as- 
sume values noticeably difi'erent from those of the FG stars. 

The analysis of the relative positions of the turn-off, 
the red giant bump and the HB could not be performed on 
other blue-HB clusters, due to the lack of a consistent pho- 
tometry for these features. Clearly a possibility is that in all 
the clusters of this group most of the first generation stars 
have been lost. There are contradictory spectroscopic indi- 
cations in favour or against this hypothesis. In the cluster 
N GC 6752, out of nine d warfs and nine subgiants analyzed 
bv lCarretta et al.l (|2005l ). only one (subgiant) has a normal 
nitrogen content, all the others are subst antially N-enriched 
([N/F e]~ 1 - 1.7. The recent new data bv lYong fc Grundahll 
(2008), on the contrary, contain also nitrogen normal stars. 
This cluster could be similar to M13, also in having a very 
helium rich population producing the extreme HB stars. 



6.2 NGC 6397, M53: apparently SSPs, but 
possibly SG-SSP? 

Notice that NGC 6397 has [Fe/H]=-2 (JGratton et al.ll200ll ) 
and so is not an intermediate metallicity cluster. It has a 
short blue HB, lacking extreme HB and blue hook stars, and 
its HR diagram has always been regarded as a perfect exam- 
ple of SSP, especially following the excee dingly refined HST 
pro per motion se lected observations by iKing et al.l (11998! ) 
and [Richer et al.l (|2006il . Nevertheles s, only three subgiant s 
out of 14 stars are nitrogen normal (|Carretta et al.l 120051 1 . 
leading us to suspect that the material from which these 
stars formed is CNO processed and thus of SG. Thi s occur- 
rence had already been noticed in iBonifacio et al.l ([2003), 
with reference to the paradox that nitro gen rich stars 
had a lmost-normal lithium content (see also IPasguini et al.l 
120081 ). 

The above considerations lead naturally to the ques- 
tion: the simple population GC - one chemical composition, 
one age, of which one used to speculate until very recently - 
does it exist? We are inclined to give a negative answer. Let 
us consider the observations bv lLi fc BursteinI (|2003l ). who 
took integrated spectra of eight galactic GCs, ranging in 
metallicity from [Fe/H] < - 2 to [Fe/H]~ -0.8, and show- 
ing a variety of HB morphologies (the clusters are: M 15, 
M 92, M 53, M 2, M 3, M 13, M 5, M 71). All these clus- 
ter show a substantial N-enhancement with respect to field 
stars of the same metallicity. Since we are dealing with inte- 
grated spectra, the result can not be directly interpreted in 
terms of percentage of second generation stars, that would 
be composed by nitrogen rich (CN or CNO cycled) matter. 
Nevertheless, the N- excess with respect to the field is a clear 
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indication that at least a certain amount of stellar matter is 
not the original one, in all the eight GCs quoted above. 

The case of M 53 presents some interesting features. 
It is a very metal poor ([Fe/H] ~ —2, Zinn 1985, Harris 
2003), massive cluster (My = —8.70 mag). Out of a total of 
307 observed HB members (|Rev et aLlliggsl ). 257 are located 
blueward of, 35 within, and 12 redward of the RR Lyr insta- 
bility strip, the bluest stars being located at (B-V) ~ -0.08. 
Three stars appear separated at bluer colours. So, the HB 
has a very short extension in colour, with most of its mem- 
bers concentrated at (B-V) ~ 0.05. This star distribution is 
unique among very metal poor GCs - see, f.e.. Table 13 in 
IWalkerl l| 19941 ). In fact: in M 15, a massive and high central 
density cluster, the 155 HB members are distributed from 
redward of the RR Lyrs to very blue colours and low visual 
luminosity, reaching the TO magnitude; M 68, a relatively 
small and not very concentrated cluster, exhibits a short HB 
from the red region to about (B-V) ~ -0.1 mag, with almost 
equal numbers of blue on the one side, and red and variable 
stars, on the other. 

The extraordinary star concentration in a small colour 
interval suggests that the largest part of the HB popula- 
tion in M 53 can be obtained assuming a dispersion in the 
mass loss of ~ 0.01 Mq; a more precise photometry could 
give more stringent limits to the dispersion. On this ba- 
sis, M 53 would appear a good candidate for a "normal", 
first-generation-only cluster, but for the N-enhancement 
and the the (mild) intrinsic spread in CN bandstrength 
IjMartell et al.ll2008r ). These chemical properties suggest that 



we are likely dealing with a one-generation cluster, but one 
in which the star generation we observe is the second one, 
and not the first. 

Of course a final answer to the initial question will have 
to wait the investigation of the whole body of Galactic GCs, 
but since now secondary episodes of star formation appear 
widespread and crucial for building-up the clusters them- 
selves. 



7 CONCLUSIONS: HOW DID THE GCS 
FORM? 

In all GCs examined in this work, a large fraction of the stel- 
lar population takes origin from secondary star formation 
episodes. Notice that we have examined only a fraction of 
the clusters with HB morphology or RR Lyr period distribu- 
tions similar to those described here, so that we can suggest 
that the results of this work probably hold for a larger pop- 
ulation of Galactic GCs. While the most massive clusters 
have extreme helium enhancements, also moderately mas- 
sive clusters show a considerable degree of helium variation. 
We reached our goals by examining in detail GCs that 
have unexplained features in their HBs, and extending the 
results to clusters with similar features. The HB morphology 
is one of the important features: clusters having a bimodal 
or multimodal HB are most easily interpreted by the coexis- 
tence of multiple generations with different helium content. 
Also a unique SG, whose stars had different degrees of mix- 
ing with pristine matter and thus ended up with different 
helium contents, is a possible solution. In any case, the ex- 
treme helium rich populations (in to Cen and NGC 2808 at 



least) are neatly separated from the other MS stars so that 
they should have a well defined independent origin. 

Further, we used the period distribution of RR Lyrs in 
several clusters in order to reject the hypothesis of a unique 
Y value with a relatively large spread of mass loss on the 
RGB, that has been the standard way of interpreting the 
whole HB colour extension, but is inconsistent with most 
period distributions. 

We re-examined in detail the HB distribution in 
NGC 2808, and obtained a helium di stribution co n sisten t 
with the main sequence recent data bv lPiotto et al.l (|2007r ). 
Also the UV data of this cluster find a good interpretation 
in terms of population with varying helium content, if we 
make the further hypothesis of deep mixing to understand 
the location of the blue hook stars. We also show that sim- 
ulation of M 13 UV data is well explained with populations 
having difi'erent helium. 

After this analysis, then, we must face the problem that 
the SG formation is not a peculiarity of a few very mas- 
sive clusters, but must be the normal way in which a GC 
is formed in our Galaxy. It is almost obvious, and has often 
been discussed in the literature, that the ejecta of a unique 
first stellar generation with a normal initial mass function 
(IMF) can not produce enough mass to give origin to such a 
large fraction of second generation stars (see, e.g., the case 
made by Bekki and Norris (2004), for the blue main se- 
quence of ojCen). The only solution to the IMF problem 
is that the starting initial mass from which the first gener- 
ation is bom was much larger than today first generation 
remnant mass (at least a factor 10 to 20 larger), so that 
the processed ejecta of the first generation provide enough 
mass to build up the second one. There are two possible 
ways of producing this result: the first possibility is that 
all these GCs formed within a dwarf g alaxy environment 
(|Bekki fc Norrisll2006l : ISekki et all 120071 ). There, GCs may 
be formed by mixing of pristine gas with the winds of the 
very numerous massive AGB stars evolving in the field of the 
dwarf galaxy, and later on the dwarf galaxy is dynamically 
de stroyed. A se cond possibility has been recently suggested 
bv lD'Ercole et a l. (2008). They show that the SG stars are 
preferentially born in the inner core of a FG cluster, where 
a cooling flow collects the gas lost by the FG stars. The 
massive stars that explode as SN II were preferentially con- 
centrated in the cluster core. After the mass loss due to the 
supernovae type II explosions, the cluster expands, and be- 
gins losing the stars — mainly of FG — going out of the tidal 
radius. Thus the cluster may be destroyed, unless the gas 
lost by the most massive AGB stars begins collecting in the 
core and forms the SG, that initially does not take part in 
the cluster expansion. The study of the cluster dynamical 
evolution, followed by means of N-body simulations, shows 
that the cluster preferentially loses FG stars; these simula- 
tions show that high SG/FG number ratio can be achieved 
and SG-dominated clusters may survive. 
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Figure Al. The HB data by Bedin et al. (2000) (open triangles) 
and their simulation superimposed (full triangles). The histogram 
of the blue HB data is shown as a full histogram on the right, 
and the simulated histogram is dot-dashed. The lower luminosity 
clumps EBT2 and EBT3 are both obtained with a unique value 
of Y=0.385, as explained in the text. 
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Figure A2. The full histogram represents the number vs. he- 
lium content distribution of the NGC 2808 HB of Fig. Al. The 
dashed histogram represents the corresponding distribution for 
M 13. Notice that the M 13 distribution completely lacks stars 
with normal Y=0.24. 



APPENDIX A: A NEW ANALYSIS OF 
NGC 2808 OPTICAL AND ULTRAVIOLET 
DATA 



Al The V vs. B-V data 



We a nalyzed the [B edin et al. ' (2000*) data of NGC 2808 
twice l|D'Antona & C aloi 2004; D'An tona et al.l 20051. but 
always assuming that the mass lost on the RGB has a slight 
but positive dependence on the helium content of the sam- 
ple. Having now shown that this is not the case (Sect. [2} 
it is reasonable to make another analysis, and derive the 
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MS Simulation Z=10-3 
17-0.012 age=llGyr 




MS Simulation Z=10-3 
17=0.015 age=llGyr 




MS Simulation Z-lO-a 
17=0.017 age-UGyr 
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Figure A3. The figures show simulation of the main sequence colour distribution corresponding to the N(Y) distribution of Fig. lA2l in 
the ACS colour F475-F814 of Piotto et al. 2007. The simulation is shown for the magnitude interval 5.0^ Mgi^ ^6.0. The three panels 
correspond to different assumed errors on the colour (0.012, 0.015 and 0.017mag). 
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^21B '^ 555 oiu ^^u 

Figure A4. We show the ZAHBs of Z=0.001 with different helium content in the absolute UV magnitude F218 versus F218-F555 colour. 
The Y=0.24 ZAHB line is used to define the distance modulus and reddening of the dataset by Castellani et al. 2005. We impose that 
the red clump data (on the right bottom part of the left figure) lie symmetrically on this ZAHB. The rising lines on the left are ZAHBs 
of models with the same core mass of the Y=0.40 sequence, but having larger envelope helium abundance. From left to right, we have 
Y=0.8, 0.7, 0.6, 0.5 and 0.45. On the right, we show the simulation of Fig. Al. 



N(Y) distribution that fits the HB, obtained by assuming 
that both 5Mo and a do not depend on Y. For the simula- 
tion, w e cl osely follow the proced ure bv lD'Antona fc Caloil 
l|2004 ) and lD'Antona et al.l (|2005l '). apart from the analysis 
of the blue hook sta rs (the clump EBT3 in the definition 
by iBedin et a l. 2004) for which we make the further assum- 
tions described below. The problem is that we can not adopt 
a different helium content for each of the clumps EBT2 and 
EBT3, because, using the N(Y) distribution that reproduces 
the HB, we have to reproduce also t he colour distribution in 
the MS data by iPiotto et al.l (|2007l ). For this aim, we con- 
vert our main sequence stellar models into the ACS bands 
F4 75 and F814 by m eans of the transformations provided 
by iBedin et al.l l|2005l ). and simulate the MS colour distri- 
bution. We can not make a detailed comparison with the 



data, as they are not available to us, but simply show the 
histogram of number versus the colour difference (in F475- 
F814) from the reference main sequence of Y=0.24, in the 
magnitude interval 5 ^ M^ii ^ 6. Compa rison with the his- 
tograms in Figure 3 in lPiotto et al.l (|2007l ) shows that there 
is a fair reproduction of the observations, if the colour error 
is taken to be 0.017mag (indicated by a in the top of Fig. 
[1. 



Let us discuss in detail the assumptions made to fit the 
EBT2 and EBT3 clumps. As the blue MS is well separated 
from the other stars, all the stars in this blue MS share the 
same helium abundance. We then take a unique high helium 
abundance for all the st ars in the clumps EBT 2 and EBT3 
(as we assumed also in iD'Antona et al.ll2005h . Attributing 
to the cluster an age of llGyr, the average mass loss rate 
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necessary to fit the red clump is 0.18 Mq. Due to tiie spread 
in mass loss (cr^O.OOS Mq) assumed in order to fit the width 
of the red clump, and due to the choice of a very high helium 
in order to reproduce the blue MS, the remnant mass is in a 
very strict range close to the helium flash mass (M^O.49 Mq, 
to be compared with Mc=0.4676Mq, for Y=0.40). We wiU 



fix a value M„ 



above which we put the star on the 



track corresponding to its mass and helium content. Below 



M. 



we assume that the star has suffered very deep 



mixing, and that its helium surface abundance has increased. 
We parametrize the resulting surface helium abundance be- 
tween fixed values, and distribute random the stars along 
the corresponding tracks. By this hypothesis, we are able to 
reproduce well the gap between EBT2 and EBT3 without 
invoking a helium discontinuity, as shown in Fig. lAll We 
summarize the parameters chosen for this fit: age of llGyr, 
AMo = O.I8M0, a = 0.008 
N(Y) is shown in Fig. lA2l 



M^i^i„3 =0.487 M0, and the 



A2 The HST UV data of NGC 2808 

The optical bands are certainly not the best to describe 
the very hot HB stars in NGC 2808. The HST data have 
shown that the se objects are the most luminous on es in 
the UV bands (JBrown et al.l bOOll l. iLee et all l|2005l ') first 
attempted to fit the data with HB model s having varying 
helium content. Here we use the data by ICastellani et al.l 
l|2006t ) in the plane F218 versus the colour F218-F555. Our 
models have b een transformed into these bands by using the 
lOriglia fc Leitherer ( 2000i ) colour transformations. Fig. IBll 
shows on the left the data and the ZAHBs. The choice of 
the distance modulus and reddening in the colour F218-F555 
have been made in order to fit the red clump data on the 
Y=0.24 ZAHB. This choice leads to a good superposition 
of the hig hest luminosity sta rs (corresponding to the EBTl 
clump bv iBedin et~aLl 120041 ') above the ZAHB of Y=0.32. 
The rising short lines on the left are ZAHBs of models with 
the same core mass of the Y=0.40 sequence, but having 
larger envelope helium abundance. From left to right, we 
have Y=0.8, 0.7, 0.6, 0.5 and 0.45. The simulation shown in 
Fig lAll for the B and V colours is shown in the right part 
of Fig. El in these UV bands. The fit of red clump, EBTl 
and EBT2 is very good, showing once for all that there is a 
helium difference between the red clump, EBTl and EBT2. 
The gap between EBT2 and EBT3 results well simulated by 
assuming that a fraction of the stars born with Y=0.40 suf- 
fers deep mixing, increasing its Y to the range Y=0.7-0.8. 
Nevertheless, the simulation fails to reproduce the lowest lu- 
minosity EBT3 stars. Reasons for this result are several: first 
of all, these stars have atmospheres which are not only very 
helium rich as we assumed, but also carbon rich, and the 
Teff-colour transformations can not take this into account. 
Second, it is well possible that these stars ignite helium with 
a late fiash at a core mass smaller than we assumed, and thus 
have a smaller intrinsic luminosity. 



APPENDIX B: A POSSIBLE FIT FOR THE HB 
OF M 13 




-1 

F336-F555 

Figure Bl. M 13 HB data by Ferraro et al.(1998) compared 
with our ZAHBs in the plane F555 vs. F336-F555. The distance 
modulus is chosen so that the Y=0.28 ZAHB (dotted line) fits 
the upper luminosity clump. The other lines are the ZAHBs for 
Y=0.24 (full line), Y=0.32 (dashed) and Y=0.40 (dash dotted). 
On the left we show the histogram of HB counts as a function of 
F555 magnitude. 




Figure B2. Superimposed to the data, we show the sim- 
ulation done with the same values, <5M=O.18M0 and age 
11 Gyr, chosen for the NGC 2808 simulation, a is 0.01 Mq and 



M„ 



:.479 Mq 



ICaloi fc D'Antonal (|2005l ) provocatively proposed that the 
blue-HB clusters, and in particular M 13, have completely 



lost their FG (see Sect. 6.1). As the metallicity of M 13 
and NGC 2808 are close, and the blue HB of the two clus- 
ters are morphologically similar, we try to fit the HB of 
this cluster by imposing a simulation in which the age and 
average mass loss are the same as in the NGC 2808 fit, 
but the red clump population, at Y= 0.24, is totally elimi- 
nated. FiglBll shows the HST data bv lFerraro et all (Il998l ) 
in the plane F555 vs. F336-F555. We superimpose our ZA- 
HBs as in Fig. IBll assuming a visual distance modulus of 
14.6mag and zero reddening. The modulus is chosen so that 
the Y=0.28 ZAHB coincides with the upper clump, in agree- 
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ment with the resuh by ICaloi fc D'Antonal (120051 ) that at- 
tributed Y=0.28 to the dominant cluster population. We 
see that, with this choice, the ZAHB of the middle HB is 
again consistent with Y~0.32. This difficulty of fitting the 
HE with a unique ZA HB had already been pointed out by 
iGrundahl et al.1 (|l998l ) in their analysis of the HR diagram 
of M 13 in the Stromgren colours, and signalled also for the 
clusters NGC 288 and NGC 6752. These authors indeed at- 
tributed this unexpected feature to the presence in M 13 
of two distinct HB pop ulations, one of which had undegone 
deep mixing, following ISweigartI (|l997l ). In our interpreta- 
tion, these stars have a higher helium content already start- 
ing from their formation. 

Fig. IB2I shows the HB simulation superimposed to the 
M 13 data. Apart from the lack of the Y=0.24 part, the 
distribution N(Y) is different indeed from that obtained for 
NGC 2808, as we show in Fig.|X2l but the lowest part of the 
HB can be interpreted again as a very high helium popula- 
tion (Y=0.38). Notice however, that M 13 does not contain 
the large population of blue hook stars present in NGC 2808, 
and this difference remains to be explained. 



